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Abstract 

Particle track reconstruction capabilities of the silicon tracking detector system 
have been studied. As the multiple Coulomb scattering (MCS) induces unavoidable 
uncertainties on the coordinate measurement, the corresponding error estimates and 
the associated correlations have been used to find the best track fit parameters and 
their errors. Finally it permits to find the proper particle characteristics, as vertex 
position and resolution, flight direction and the error. 
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1 Introduction 

Design and preparation of any components of a detector system must take care of charac- 
teristics and detection performances (efficiency, acceptance, position or energy resolution) 
necessary for a specific process study. The silicon tracking detector system, in our case, 
must furnish the best information about the coordinate track intercept of the incident par- 
ticle on every detector layer. Otherwise, for a good track reconstruction, it is necessary 
to find the best estimate of the track parameters in a specific particle measurement. This 
is the reason we studied, the particle transport in a silicon tracking system and estimated 
the multiple Coulomb scattering (MCS) perturbation in particle track measurements. 

This study furnishes information on the position resolution we could get from every 
detector layer, and also its dependence on some physical and geometrical parameters. In 
the detector system design these parameters could be choosed in an optimal manner, aim- 
ing the best track reconstruction possibilities or the best estimate of the track parameters 
(vertex position, curve radius) directly connected to physical quantities. 

The simplified tracking system we used, consists of a 5 layers silicon microstrip detec- 
tor, 300/im thick, interspaced by 1.5 cm, and the first detector layer located at 130 cm 
from the interaction point (see Fig.l). 

2 The track particle position uncertainties due to 
multiple Coulomb scattering. 

When a charged particle is traversing the detector elements of a tracking system, it un- 
dergoes small deviations of the track, due to MCS. The effect is usually described by the 
theory of Moliere (see for example |T[). It shows that, by traversing detector's material, 
thickness s, the particle undergoes successive small-angle deflections, symmetrically dis- 
tributed about the incident direction. Applying the central limit theorem of statistics to a 
large number of independent scattering events, the Moliere distribution of the scattering 
angle can be approximated by a Gaussian one [H. It is sufficient for many applications 
to use Gaussian approximation for the central 98% of the plan projected angular distri- 
bution. The width of this distribution is the root mean square of the scattering angle 



where p, (3c and z c are the momentum, velocity and charge number of the incident particle, 
and Xl is the radiation length of the scattering medium. That is, the plane projected 
angle 9 p i aneyX or 9 p i ane ^ y of the deflection angle 9, onto the xOz and yOz planes, where 
the x and y axes are orthogonal to the Oz direction of motion, shows an approximately 
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Deflections into Q p i an e,x and Q v iane, y are independent and identically distributed, and 6 2 space = 

plane, x ' plane,y 

The angular distribution is translated to a coordinate distribution by particles fly onto 
every detector layer. The more intersected layers the larger distribution width is. The 
coordinate distribution is defined by statistical spread due to MCS, and depends on the 
number and position of the intersected detector layer elements. It has the same form as 
angular distribution 
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with the mean square deviation (distribution width) as the squares sum of the (i — 1) 
preceding distribution widths projected onto i-th detector layer 
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For an oblique incidence (8 ^ 0) the effective path length in the silicon detector is larger, 
and the same is the position distribution width on the next detector layers. Nevertheless, 
in this work we will consider only the minimal width (f4|), as to emphasize the precision 
limit in particle position measurement with a given silicon tracking system. 



3 The Monte-Carlo particle scattering description. 

The change of track parameters is usually parametrized by two mutually orthogonal, 
uncorrelated scattering angles 59 p i anejX and 56 p i ane ^ y or (59 space and 5(p) which leads to a 
corresponding displacement Sx p i ane and Sy p i ane in the position. 

Following the stochastic nature of the MCS, we use the Monte-Carlo study by generat- 
ing the joint (5x p [ ane , S9 p i anejX ) distribution with independent Gaussian random variables 
(wi,w 2 ) § 

WxsQq w 2 s9 

OXplane ~ 2 
&Qplane,x = W 2 0q (5) 

The same have been used for the joint {5y p i anei 56 p i ane ^ y ) distribution in yOz plane. 
Finally we constructed the incidence points distribution (coordinate distribution) on the 
detector layer no. 2, 3, 4 and 5. (see Fig. 2 and Table [1|), for an incident 500 MeV/c muon. 
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In Monte-Carlo simulated particle transport the particle position uncertainty on every 
detector layer has been measured as the variance or the mean square deviation of the 
scattered track incidence points x m about the unscattered one xq 

N f„ „ \2 
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and similar for <jy variance, for every detector layer. The x n (and y n ) are the coordi- 
nates of the scattered track incidence points on detector layer, and N is the total number 
of generated events. 

The distribution widths a Xi of the muon track x-coordinate points on detector layer 
i= 2, 3, 4, 5 and for momentum values from 500 to 7500 MeV/c, are presented in Table [1]. 
There are also the analytical estimation of the same widths of the MCS incidence point 
distribution on each of the i-detector layer, according (|J) 

< = ^f-V )' (7) 

where / - the distance between detector layers (/ = 1.5 cm) and #o - the plane r.m.s. 
scattering angle ([!]). 

This Monte-Carlo particle transport description within silicon detector system will be 
used as a coordinate data generator for the track reconstruction procedure. 



4 Position error correlations. 

MCS produces errors correlated from one layer to the next. Clearly a scattering in layer 1 
produces correlated position errors in layer 2, 3 and so on (see Fig. 3). The proper error 
matrix is non-diagonal |J, and it must be find out. 

Let's denote Sxi the track deviation x-coordinate point on the i-th layer, with respect 
to the initial incident direction on the detector system, then (see Fig.3) 

6xi = 
6x 2 = 6 l 

6x 3 = 5 31 + 5 2 (8) 

6X4 = 641 + 642 + 63 

6x 5 = 6 51 + ($52 + <fes + &t 

where the individual contributions due to preceding scatterings are 
631 = 61 
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Because they are statistical variables it is necessary to find their mean value, and to 
express it by the independent scattering deviations 5k, unaffected by earlier scatterings 



Si 



<si> 



(10) 



We try to express the 5xi deviations by the independent 5k ones (|10"D, for which 

< 5 k >= 



r- V~~V r- Z{ Zk 

5xi = 2^ 5 k 

k=i z k+i ~ z k 



Now the position error (covariance) matrix (Vij), defined || as the statistical mean of 
the pair deviation products < 5xi5xj > for all possible detection layers is 



Vij =< 5xi5xj >- 



9l [{zi - z l )(z j - zi) + (zi - z 2 ){zj - z 2 ) + ... + {zi - Zi-i)(zj - Zi-i)] (12) 



for i<j= 1,2, ...,n 



The {Vij) matrix is symmetric. The error correlation matrix is immediately 
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The uncorrelated position errors in the coordinate reading Co, have to be added in squares 
(<7q) into the diagonal terms of the error matrix V. 

For example, for a 500 MeV/c MUON track detection (6* = 1.23mrad) by our system 
configuration (z± = 130cm, z 2 = 131.5cm, z% = 133cm, z^ = 134.5cm, z$ = 136cm and 
o"o = 10/im) the matrix elements Vij (in square microns) and pij are 
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( 100. 
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0. 

V o. 
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1021.9 
1362.5 
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1803.2 
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4868.9 
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0. \ 
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10319. J 
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.63898 
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.96114 
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(14) 



/ 1 

1 .76430 .69768 

p = .76430 1 .91970 .86865 (15) 

.69768 .91970 1 

v .63898 .86865 .96114 

As long as the Vij matrix elements depend both on kinematical characteristics of the de- 
tected particles and on the tracking detector system configuration, the matrix elements 
are independent on particle characteristics, and is defined only by system configuration. 

In the following we will use these matrices in the track reconstruction by a least squares 
fit procedure. 
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5 Track reconstruction parameters and the errors 



In the absence of the magnetic field, the unscattered track is a straight line. The inde- 
pendent description of x and y MCS data permits a separate least squares fit to these 
data by a linear relationship || |7j 
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With the coordinate and error data (xj ± <t x J, (?/« ± a Vi ), z iy along with the corresponding 
correlation matrix (Vij = PijCr Xt a x .) as input data, it is possible to express the \ 2 m 
matrix form, for every coordinate data set. For x-data set it will be 
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By solving the linear system relative to A x we get the fit parameters 

A x = {H T V- 1 H)- X {H T V- 1 X) 
and the errors of these parameters 

E Ax =< 8AJAZ >= (H^-'H)- 1 



(19) 



(20) 



The same procedure we applied to y coordinate, as to find the fit parameters and errors; 
Zi has no any uncertainty, they are detector position coordinates. 

Let's take an example. Using the particle transport simulation data within silicon 
tracking system, along with the variance matrix data V (|I4D , we have been obtained the 
fit parameters and their errors (x ± <j xo , a x ± a ax ; y ± a yoi a y ± cr ay ). The results are 
presented in Table 0. 

The found track fit parameters, could be connected with some particle characteris- 
tic quantities as xy-vertex position (xq,?/o) and the xy-vertex uncertainty (<J x v y ertex ^ = 
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\J a x Q + a y ) an d the particle production direction {9 ± a e ,(p ± a v ). This last one is con- 
nected with the slope parameters a x and a y . See Table 0. 
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and the corresponding errors 
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aL = I — I ^ + I — ^- I cr* (22) 



't<70 

In Table [| there are results on the reconstructed and ip direction values in comparison 
with the generated ones. 

The xy-vertex position unc ertainty, due to MCS in the detector material, have been 
also calculated as a^ y ertex ^ = \J<?1 + a y an d depends, of course, on particle momentum. 
If we want to reduce this error it is necessary to bring closer the detector tracking system 
relative to interaction point. The xy- vertex position uncertainty (resolution) dependence 
on particle (muon) momentum and z\ distance to first detector layer is shown in Fig.4. The 
result is useful in detector system design to find an optimal configuration in preparing the 
experimental work. Also, to have a choice for vertex resolution as a compromise between 
the best possible values in the proximity of the interaction point and the worse ones far 
from this point at the radiation harmless distance. The xy-vertex position resolution due 
to MCS errors combined with the intrinsic detector coordinate uncertainties, could not 
be better than the values shown in Table |j. Nevertheless, from Table f| and Fig.4 we see 
also that it is possible to have a better vertex position resolution if the detector tracking 
system is placed to a smaller distance from the interaction point. 



6 Conclusions 

We have calculated the uncertainties in particle characteristic parameters as vertex po- 
sition and the particle flight direction, by a linear fitting procedure, using position error 
correlation due to MCS in succesive detector layers of a silicon tracking system. 

The most important result shows these uncertainties are independent on individual 
track coordinates (see for example the errors for x ,y , and a x ,a y ) and they depend exclu- 
sively (see (|2CiD and (|^)) on the detector system configuration (geometry). The results 
are useful both in particle tracking system design and in experimental data analysis for 



pattern track recognition. The former one could use our analytical relations for parti- 
cle parameter error estimation relative to diferent tracking system configuration. The 
later one could use our data fitting procedure, including system configuration dependence 
on the coordinate covariance matrix, in particle track reconstruction and the associated 
parameter determination. 
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Table V. MUON x-distribution width on the detector layers of a silicon tracking system 
due to multiple Coulomb scattering. 

These are the correlated uncertainties <7j, used in track fitting procedure. 

MONTE-CARLO SIMULATED DATA / ANALYTICAL CALCULATED DATA (g) 



p 






a 2 








0"3 






04 
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(MeV/c) 




(fim) 










(/im) 




(/im) 


500 


18 


88 


1 


1 


3.46 


41 


75 


/ 


41.27 


69 


74 


/ 


69.06 


101 


94 


/ 101.09 


1000 


9 


20 


I 


9 


08 


20 


60 


/ 


20.30 


34 


43 


/ 


33.97 


50 


34 


/ 49.73 


1500 


6 


10 


/ 


6 


03 


13 


55 


/ 


13.49 


22 


73 


/ 


22.58 


33 


42 


/ 33.05 


2000 


4 


56 


I 


4 


52 


10 


19 


/ 


10.11 


17 


05 


/ 


16.91 


25 


07 


/ 24.76 


2500 


3 


71 


/ 


3 


61 


8 


18 


/ 


8.08 


13 


65 


/ 


13.53 


19 


96 


/ 19.80 


3000 


3 


08 


/ 


3 


01 


6 


84 


/ 


6.73 


11 


36 


/ 


11.27 


16 


57 


/ 16.49 


3500 


2 


62 


/ 


2 


58 


5 


84 


/ 


5.77 


9 


72 


/ 


9.66 


14 


18 


/ 14.14 


4000 


2 


28 


1 


2 


26 


5 


09 


/ 


5.05 


8 


54 


/ 


8.45 


12 


52 


/ 12.36 


4500 


2 


03 


1 


2 


01 


4 


50 


/ 


4.49 


7 


48 


/ 


7.51 


10 


90 


/ 10.99 


5000 


1 


82 


1 


1 


81 
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05 
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4.04 


6 


79 


/ 


6.76 


9 


96 


/ 9.89 


5500 


1 


66 
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1 


64 


3 


68 
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3.67 


6 


12 


/ 


6.14 


8 


94 


/ 8.99 


6000 


1 


53 


1 


1 


51 


3 


39 


/ 


3.37 


5 


67 


/ 


5.63 


8 


31 


/ 8.24 


6500 


1 


39 


/ 


1 


39 


3 


12 


/ 


3.11 


5 


22 


/ 


5.20 


7 


61 


/ 7.61 


7000 


1 


31 


/ 


1 


29 


2 


90 


/ 


2.88 


4 


85 


/ 


4.83 


7 


06 


/ 7.07 


7500 


1 


21 


I 


1 


20 


2 


69 


/ 


2.69 


4 


50 


/ 


4.50 


6 


58 


/ 6.59 



Table 2: Particle transport simulation data within silicon tracking system. 



MUON 500 MeV/c, x = 0, y = 0, z = 0, 9 = 15°, <p = 30 
Silicon layer width = 300/xm, z\ = 130cm, dz = 1.5cm 
Mean scattering angle (plan projected) 9 = 1.23mrad 



det 


x(cm) 


a ((j,m) 




y(cm) 


cro(fjim) 


crMCs(l^m) 


z{cm) 


1 


30.166605 


10 


0.0 


17.416698 


10 


0.0 


130.0 


2 


30.517353 


10 


18.5 


17.617701 


10 


18.5 


131.5 


3 


30.868230 


10 


41.3 


17.816364 


10 


41.3 


133.0 


4 


31.218312 


10 


69.1 


18.013349 


10 


69.1 


134.5 


5 


31.569294 


10 


101.1 


18.207651 


10 


101.1 


136.0 



Linear track reconstructed (parameter) data 

x =-2316.4 ± 1828.3/im a x =0.23383 ± 0.001404 
y = 618.69 ± 1828.3^m 0^=0.13350 ± 0.001404 

a (vertex) = = 2 585.65fim 

= 15.070° ±0.096° 
up = 29.723° ± 0.298° 



Table 3: Track reconstructed direction by least squares fit of the 500 MeV/c MUON 
Monte-Carlo simulated transport data within silicon tracking sistem. 



Generated 


Reconstructed 


0(deg) 


(p(deg) 


9(deg) 


(p(deg) 


5 


270 


5.07 ±0.81 


269.20 ±0.91 


10 


240 


10.16 ±0.20 


239.52 ±0.45 


15 


210 


14.87 ±0.10 


210.20 ±0.30 


20 


180 


19.85 ±0.07 


179.79 ±0.22 


25 


150 


25.02 ±0.08 


149.89 ±0.17 


30 


120 


30.05 ±0.16 


119.85 ±0.14 


35 


90 


35.04 ±0.20 


89.98 ±0.11 


40 


60 


39.91 ±0.13 


59.94 ±0.10 


45 


30 


44.99 ± 0.05 


30.04 ±0.08 



Table 4: xy-vertex position uncertainty due to multiple Coulomb scattering by track 
reconstruction procedure. 



Table shows the xy-vertex uncertainty (/im) obtained by a linear track reconstruction 
with a 5 detector layers system, interspaced by dz—1.5 cm, for some z\ target to first 
detector layer distance (cm) and some MUON momentum (MeV/c) values. 
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203.414 


600 


261 


997 


307 


1394.383 


1791.468 


2188.558 


2585.650 


1000 


116.415 


336 


970 


557 


795 


778 


662 


999 


541 


1220.428 


1441.317 


1500 


86.285 


245 


645 


405 


317 


565 


036 


724 


772 


884.515 


1044.262 


2000 


70.975 


198 


945 


327 


254 


455 


615 


583 


994 


712.382 


840.774 


2500 


61.993 


171 


278 


280 


916 


390 


610 


500 


322 


610.044 


719.770 


3000 


56.277 


153 


481 


251 


041 


348 


659 


446 


296 


543.943 


641.594 


3500 


52.429 


141 


375 


230 


675 


320 


033 


409 


410 


498.797 


588.189 


4000 


49.726 


132 


790 


216 


202 


299 


672 


383 


161 


466.660 


550.164 


4500 


47.761 


126 


498 


205 


575 


284 


709 


363 


863 


443.026 


522.194 


5000 


46.291 


121 


760 


197 


561 


273 


417 


349 


292 


425.177 


501.067 


5500 


45.166 


118 


111 


191 


379 


264 


701 


338 


043 


411.394 


484.749 


6000 


44.287 


115 


246 


186 


519 


257 


846 


329 


191 


400.546 


471.905 


6500 


43.588 


112 


957 


182 


634 


252 


363 


322 


110 


391.867 


461.628 


7000 


43.024 


111 


103 


179 


482 


247 


914 


316 


364 


384.822 


453.285 


7500 


42.562 


109 


581 


176 


894 


244 


258 


311 


640 


379.031 


446.426 



Figure captions 



Fig. 1 The track particle position uncertainties due to multiple Coulomb scattering. The 
uncertainty is defined as the coordinate distribution width <7j on every detector 
layer i and results as the squares sum of all preceding layer scattering contributions 
projected onto the i-th layer, a\ = ^[(^ ~ z i) 2 + ( z i ~ z i) 2 + ••• + ( z i ~ z i-i) 2 ] 

Fig. 2 The plan projected x-coordinate point distributions of the scattered 500 MeV/c 
MUONS, incident on detector layer No. 2, 3, 4, and 5 obtained by Monte-Carlo 
particle transport simulation. 

Fig. 3 Multiple Coulomb scattering error correlations. 

Fig. 4 The xy-vertex position uncertainty (reconstructed vertex resolution deterioration) 
due to multiple Coulomb scattering as a function of MUON momentum and distance 
from the interaction point to first detector layer (zi). 



